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O 

f"*^ ' We consider the Coulomb blockade on a superconductive quantum dot strongly coupled to a lead 

£\j | through a tunnelling barrier and/or normal diffusive metal. Andreev transport of the correlated pairs 

leads to quantum fluctuations of the charge on the dot. These fluctuations result in exponential 
renormalization of the effective charging energy. We employ two complimentary ways to approach 
the problem, leading to the coinciding results: the instanton and the functional RG treatment of the 
non-linear sigma model. We also derive the charging energy renormalization in terms of arbitrary 
OO ' transmission matrix of the multi-channel interface. 
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PACS numbers: 73.50.Bk, 72.10.Fk, 72.15.V, 73.23.Hk 



ryj ' I. INTRODUCTION 



Physics of interacting electronic systems in presence of disorder has been a subject of an intense study already for a 
few decades ||] ■ Various theoretical approaches have been developed for the description of both metallic and insulator 
phases. The non-linear cr-model (NLcrM) in the replica ]^J| or dynamic (Keldysh) formulation has proven to be 
the most powerful tool to deal with the weakly disordered (metallic) phase. It was shown that both versions of the 
cr-model may effectively treat the perturbation theory as well as the renormalization group (RG) formalism. Unlike 
the non-interacting models, where a whole specter of non-perturbative results is available ||, it was relatively little 
progress in the development of non-perturbative solutions of the interacting NLcrM. Our goal is to make a step in 
this direction, using the Coulomb blockade (CB) on a superconductive (SC) quantum dot as a prototypical example. 

The Coulomb blockade on a quantum dot coupled to a certain number of leads proved to be an extremely rich 
and fascinating phenomena both theoretically and experimentally (see Refs. |7],|8| for review). From the theoretical 
point of view it provides a model, where the Coulomb interactions, being spatially localized, may be treated in a 
non-perturbative way. The interactions strongly affect charge (and spin) fluctuations between the dot and the leads, 
| that manifests in the peculiar transport and thermodynamic behavior of the coupled dot-lead system. In the case 
of practically isolated dot, such that the dimensionless conductance (measured in units of Gq = e 2 /(2-Kh)) of the 
dot-lead interface is small (G <C 1), the fluctuations usually may be taken into account perturbatively. (The notable 
exceptions are provided by the vicinity of the charge degeneracy point and by the Kondo effect on the dot ||). Here 
we concentrate on the opposite scenario of the dot strongly connected to the leads (G > 1). In this case the Coulomb 
blockade is expected to be suppressed by the charge fluctuations and the overall effect of interactions on a single dot 
to be weak. It is, however, a challenging theoretical problem, lacking an obvious small parameter, to understand a 
I ■ remnants of the CB on a strongly connected dot. More importantly the weak CB may prove to be a strong phenomena 
in the granulated systems, with many (superconductive) dots incorporated in the conducting matrix. 

For the normal dot the weak CB effect was relatively well understood from the various points of view. Matveev || 
gave a complete picture of the phenomena for the case of one or two conducting channels connecting the dot and the 



leads. For the multi -channel system the CB suppression was calculated with the exponential accuracy employing RG 



technique |Py,hl|, the instanton calculus [12l-15| and bosonization jS|. The factor, suppressing the thermodynamic 



CB oscillations, was formulated by Nazarov |14| in terms of transmission coefficients, T a , of the dot-lead interface as 



]^[ Q (1 — To) 1 / 2 , where a = 1 . . . N and TV is the number of spinless channels. The remarkable feature of this result 
is that the essentially many-body phenomena may be described via the single-particle (non-interacting) scattering 
matrix only. Yet the knowledge of the interface conductance G — J2 a ^ a a l° nc is n °t sufficient to describe the CB 
oscillations. For example, in a dot coupled to a lead through the tunnelling barrier with T a <C 1, one obtains for the CB 
suppression factor exp{— Gt/2}. Another important case of a dot-lead interface is a coupling via a piece of diffusive 
metal. Employing Dorokhov statistical distribution |l(| of the transmission coefficients P(T) — Gd/{2T\/1 — T), 
where Gd is the conductance of the diffusive area, one obtains for the typical CB suppression factor exp{— tt 2 Gd/8}. 
Since the result is exponentially sensitive - the difference between e.g. G/2 and 7r 2 G/8 may be actually enormous. 

The superconductive dot in contact with normal leads is even more challenging system. We consider the low 
temperatures T <C A, where A is the SC gap, and therefore only very few (if at all) quasi-particle excitations are 
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allowed to leave or enter the superconductor. The dominant mechanism of the charge transfer through the interface 
is thus the Andreev transport of the correlated pairs. The Andreev transmission of a given channel is known to be 
T a = T£/(2 — T a ) 2 [jl7 18|. For the tunnelling barrier setup (T a <C 1) this leads to an overall Andreev conductance 
that scales like Ga ~ G\/N (pair tunnelling probability). In most cases this is a very small number. The presence 
of the diffusive normal metal adjacent to the tunnelling barrier increases the Andreev conductance to Ga ~ G t /Gd 
(or, in the case Gt > Gd, to Ga ~ Gd)- The physical reason of this phenomena is multiple attempts of Andreev 
transmission due to the coherent back-scattering on the normal impurities. The natural question is whether it is G 
or Ga (or may be none of them) that determines the amplitude of the CB oscillations. In case Ga is the relevant 
quantity one may wonder whether the coherent back-scattering enhancement should be taken into account. The 
answers are not immediately obvious, since while the SC strongly prefers the pair transport - the Coulomb energy 
of the dot makes the entrance of two charges at once energetically costly. The interactions of the dot may also 
provide a dephasing mechanism which ruins coherent back-scattering. We show that parametrically it is indeed Ga 
(including diffusive enhancement) which determines the CB amplitude. The coefficient in the exponent, however, is 
setup dependent, making the final answer sensitive to the ratio Gt/Gd and not only to Ga- The effective charging 
energy turns out to be proportional to ria(l — To) 1 / 2 , where T a are the Andreev transmissions defined above. This 
expression provides a remarkable analogy between the normal and Andreev results. 

The other interesting phenomena brought by superconductivity is the parity effect ||l9| , po| . At very small temper- 
ature the pair tunnelling is the dominant mechanism and the period (in the gate voltage) of the CB oscillations is 
twice larger than in the normal dot. For a closed dot the parity phenomena is destroyed by the entropic effects at 
moderately small temperature T* — A/ln(A/<5) -c A, where S is the mean single-particle level spacing on the dot 
pp| . The physical reason for this temperature to be much less than A is that it is enough to have a single excited 
quasiparticle to destroy the parity effect. At larger temperature the system exhibits the "normal" oscillation period. 
We show that for an open dot the transition temperature between the normal and doubled period is somewhat larger 
than for a closed dot and is given by = A/ 1ii(A/Ga<5), provided Ga& < A. 

Technically we treat the problem from the two seemingly distinct perspectives. First we look for the spatially- 
dependent instanton solution of interacting NLctM on SC dot in contact with the normal diffusive region. The finite 
action of such instanton configuration results in the exponential suppression of CB amplitude. Alternatively we can 
treat the problem employing the functional RG technique. The latter approach encodes the entire dot-lead interface 
in an (infinite) set of coefficients, whose values are subsequently renormalized by the quantum fluctuations of the 
phase on the dot. Renormalization is terminated at the cut-off energy scale where the conductance reaches unity. 
The renormalized (exponentially small) cut-off energy is the effective charging energy Ec , this dictates the amplitude 
of the CB oscillations. One of the messages of the present paper is that these two approaches lead to the identical 
results. 

The paper is organized as follows: in section [n] we describe the setup and formulate the appropriate action in the 
interacting NLcrM language. Section HJ is devoted to the instanton treatment of the problem in the real space. In 
section |V| we derive the proximity action functional and obtain CB amplitude for an arbitrary set of transmission 
eigenvalues (using instanton approach) at moderately low temperatures Ec «T« Ec ■ Then in section ^ the RG 
approach is employed to treat the CB in the same temperature range. It is demonstrated that the results coincide 
with those obtained by instanton techniques. Finally in section ^ we discuss the physical results and their possible 
experimental signatures. Appendix contains derivation of the instanton action starting from the real-time Keldysh 
functional technique. 



II. PROBLEM SETUP AND ACTION 



We consider a large diffusive (or chaotic) SC dot. The mean single-particle level spacing of the dot, S, is supposed 
to be the smallest energy scale in the problem. The SC gap, A, on the other hand, is the largest scale. The electrons 
on the dot interact via the capacitive interaction of the form 

H int =E c (N~q) 2 , (1) 

where N is the electron number operator and q is the rescaled dimensionless gate voltage potential. The charging 
energy, Ec = e 2 /(2C), is assumed to satisfy the inequalities 5 < Ec < A. 

The dot is separated from the normal diffusive metal by the tunnelling barrier with the conductance Gt, see Fig. [I]. 
The piece of the quasi ID or 2D metal having size L and conductance Gd is in turn connected to a clean bulk 3D 
lead. By the reasons which are explained later we shall assume that the Thouless energy of the diffusive region, 



2 




FIG. 1. (a) Schematic view of a SC dot connected to a bulk lead through a tunnel barrier with the conductance Gt and 
a piece of diffusive metal with the conductance Gd- The charging energy Ec controls the coupling between the dot and the 
gate, (b) A possible 2D realization. 



-Etii = D/L 2 , where D is the diffusion coefficient, is larger than the charging energy of the dot, £rh > Ec- The 
opposite limiting case requires a separate treatment and will be presented elsewhere. 

We shall be interested in those thermodynamic characteristics of the dot, which depend in an oscillatory way on 
the gate voltage potential, q. Specifically, we look for the free energy F(q) = —TlnZ(q), where Z{q) is the partition 
function. Both of these quantities depend on a particular realization of disorder in the diffusive region and inside the 
dot. We shall therefore look for the disorder averaged quantity like (F(q)). For an open dot F(q) is an oscillatory 
function, with a phase sensitive to a disorder realization. It is therefore convenient sometimes to calculate also the 
correlation function (F(q)F(q')) |15| ]. The latter carries information about the typical F(q), rather than the average 
one (that may spuriously vanish due to phase randomness). 

The disorder averaging may be performed in two ways either by introducing p replica of the system [^,^| and 
sending p — > in the end, or by dealing with the dynamical Keldysh formulation ||||. We shall employ both of these 
approaches to demonstrates what they are consistent and interchangeable. In the replica formalism, the NLuM is 
formulated in terms of the matrix field theory of the matrix field QfJ^r, t, t'), where a,b = 1,2... 2p are the replica 
indexes (one needs 2p replica to describe the correlation function (Z p (q)Z p (q'))) and i,j = 1,2 are Gorkov-Nambu 
indexes. The correlation function of the free energies is given by {F(q)F(q')) = \\m p ^p^ 2 ({Z p (q)Z p (q')) — 1). The 
Q-matrix obeys the constraint Q 2 = 1 and the fermionic anti-periodic boundary condition in both of its imaginary 
time arguments t,t' € [0, j3\. The coordinate r runs over the volume of the SC dot and the normal diffusive region. 
The matrix Q-field describes dynamics of the electrons; it is coupled to a scalar bosonic vector field $ a (r) which 
originates from the Coulomb interactions on the dot. Since we restrict ourselves by the simplest capacitive interaction 
Eq. ([!]), the field $ a {j) is space independent throughout the dot and vanishes outside the dot. As a result, only 
spatially independent (q = 0) component of the superconductive Q-matrix of the dot appears to be coupled to the 
Coulomb field Therefore we assume that the Q-matrix is spatially constant inside the dot, while it may have a 
non-trivial spatial structure inside the normal diffusive region. The effective action for our geometry contains three 
terms 

S = Sdot + St + Sd , (2) 

where Sdot and Sd are the effective actions of the isolated SC dot and normal diffusive region correspondingly, while 
St is the tunnelling action which provides the coupling between the two. We shall examine separately these three 
contributions. 



A. Action of the dot 



According to our model - the dot is the only region, where electrons interact via the Coulomb interaction Eq. (Q) . 
As a result, the dot's action contains the two coupled fields: ®h{t) and Qg b (T,T r ). The first one originates from the 
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Hubbard-Stratonovich decoupling of the Coulomb term pu , while the latter is the result of disorder averaging and 
integration of the fermionic degrees of freedom on the dot [H. Both of these fields are spatially independent, reflecting 
the fact that the Thouless energy of the dot is assumed to be large. The action takes the standard by now form: 




SdotlQs, *] = dT T# - - J Tr { - ^ + A ) Qs} ■ (3) 



Here the replica vector of gate voltages, q a , is defined in a way that q a = q for a G [l,p] while q a = q' for a e [p+1, 2p]. 
The SC order parameter A = (A a <r + — A a <r^)5 ab , is the same in all replica and is written as a matrix in the Gorkov- 
Nambu space, where a± — (<j\ ± 10-2) /2. Notice, that the gate voltage dependence of the partition function originates 
entirely from exp{i7r^ Q q a W a }, where wW a = Jdr &a(T~) = & a (f3) — ^a(O) is the zeroth Matsubara component of the 
$ Q (r) field. 

Since all the energy scales we consider are larger than S, one may evaluate the second term in the dot's action (|3|) in 
the saddle point approximation over the field Qs- In so doing, one disregards the mesoscopic conductance fluctuations 
of the dot-lead interface || . The saddle point value of the Qs field is given by the Gorkov Green function gauged by 
the phase $(r): 

Qf( T , T ') = e^^Afir, T ')e-^*b(r') _ ^ 

The Gorkov Green function, As, has the standard form, which is more familiar in the Matsubara basis (we assume 
the phase of the SC dot without the Coulomb interactions to be zero) 

Af M =SS nm (^ n *± ), (5) 



sm V n — cos ( 



where 

cos6>„ = — ; £ " = ; sm6>„ = — , = = (6) 

and e n = 7rT(2n + 1) is the fermionic Matsubara frequency. The phase rotation (|4|) preserves the fermionic anti- 
periodic boundary conditions if all W a are even integers. At small temperatures, T« A, however, one has cos0„ « 
and sin#„ w 1 and therefore the Gorkov matrix, As is (almost) off-diagonal in the Nambu space and local in time. As 
a result, only 2$(r) participate in the phase rotation, Eq. (H|). Therefore the odd-integer W a preserve the fermionic 
boundary conditions, as well. The odd-integers W a result in the doubling of the period of the F(q) function with 
respect to the normal case. This reflects the fact that the pair transfer is the dominant mechanism of the charge 
exchange. 

In fact, one has to be more careful and recall that, according to Eq. (||), one has to perform integrations over all 
Matsubara components of $ a (r) fields All non-zero Matsubara components may be eliminated by the gauge 

transformation, Eq. (^). As a result, they have no effect on the thermodynamic of an isolated dot at all (they are of 
major importance, of course, once the dot is coupled to the leads). The remaining (usual) integral over the zeroth 
Matsubara component, W a , must be performed explicitly. To this end one notices that nTW a enters the action as 
an imaginary chemical potential in the replica a. The (replicated) free energy of an isolated dot is thus a periodic 
function of each of W a with the period 2 (indeed the chemical potential always enters as exp{fj,N /T}, since the 
number operator, A, has only integer eigenvalues - the periodicity is apparent). The free energy possesses deep 
minima at even integer values of W a with the quadratic behavior in their vicinity Fdot{W) w ir 2 T 2 W 2 / (26), where 
6 = (d 2 "Fdot I '9 '/i 2 ) -1 is the mean level spacing of the dot. For a sufficiently large dot, where 8 < T, the integrals over 
W a may be performed in the saddle point approximation, which results in the even-integer quantization of W a . In a 
SC dot there are additional minima at odd-integer values of W a - Consequently the integration over W a is substituted 
by the summation over all integers. At T > the additional minima of the free energy at odd integer W a are not as 
deep as at the even integers, reflecting the fact that the addition of an odd number of electrons is possible by creating 
a quasi-particle. We shall evaluate now the action cost of the odd minima with respect to the even ones. 

For the even values, say W a = 0, one substitute the saddle point solution, Eq. (Q), into the action, Eq. (^), and 
obtains (we disregard for a moment the first term in Eq. (^)) 



S dot (W a =0) = ~ Tr{(a 3 d T + A)A s a } = ]T y/'l + l^l 2 • 

n 



(7) 
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This sum is divergent. However, it is only the difference of the action between different replica that has a physical 
significance. The later quantity is convergent as we shall see momentarily. For the odd integers, say Wb = 1, there is 
an imaginary component inT of the chemical potential in Eq. (^|) . It may be eliminated by the gauge transformation 
which converts the antiperiodic boundary conditions for the fermions into the periodic one. As a result, one arrives to 
the same expression as Eq. (H) with the fermionic Matsubara frequency e n = irT(2n + 1) substituted by the bosonic 
one u>„ — 2ixTn: 



S d ot {W b = 1) = - 



2tt 



E 



(8) 



Employing the Poisson summation formula, one obtains for the difference 

IA 



A(T) = S dot (l) - S dot (0) = 



E 



21 + 1 1 l V T 



4^2^^ _ |A|/T 
6 



(9) 



where the last equality assumes T <C A. In the opposite limit of the normal dot, T > A, one finds A(T) w ir 2 T/5, 
reflecting the fact that in a normal dot the odd-integer minima are absent. Although these minima persist up 
to T w A, their contribution is exponentially suppressed at T > T* , where T* is determined from the condition 

s dot (i) - s dot (o) « i MM 



ln|A|/<5 



« IAI . 



(10) 



As a result, there is an important temperature dependence associated with the odd-integer winding numbers at the 
scale T T* <C |A|. We shall see below that for the dot strongly coupled to the leads, the corresponding temperature 
scale is slightly different. 

We summarize now our discussion of the action of a large SC dot with the Coulomb interaction at T <gC A. The 
scalar potential in each replica obeys the boundary condition 



$„(/?) -$ a (0)=7rW a: 



(11) 



where W a is an integer winding number. The corresponding saddle point value of the Qs-matrix field is given by 
Eq. (H). The action of the dot takes the form: 



2p 
a=l 



-iirq a W a + A(T)S Wa 

mod 



2.1 



dT- 



$ 2 

^ a 

4E C 



(12) 



where A(T) is given by Eq. 



B. Tunnelling barrier and diffusive region action 

The tunnelling action couples the Qs-field on the dot with the Q(r = 0)-field at the point r = adjacent to the 
tunnelling barrier from the normal metal side. It has the standard form 

S T = Tt{Q s Q(0)} , (13) 

where Gt is the tunnelling conductance. 

The action of the normal diffusive region also has the standard form 0|| 

L 

S D = ^ JdrDTr{(VQ(r)) 2 }, (14) 
o 

where D is the diffusion constant and v the density of states (per spin) of the normal region. The total conductance of 
the normal region is given by Go = AirvD/L for the quasi ID geometry of the normal region and Gd = 87r 2 i/Z)/ ln(L/d) 
for the 2D geometry. Here L is the length (radius) of the ID (2D) region and d is the radius of the SC dot. We have 
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omitted the frequency term i/Tr{eQ} on the r.h.s. of Eq. (|lj), because of the assumption that -Exh > Eq > T. At 
the point where the diffusive region is attached to the normal bulk lead one has to impose the boundary condition 

Q(r = L)=Ajv, (15) 

where A^(n, m) = a3d ab 5 nm sign(e n ) is the appropriate Q-matrix of the normal bulk lead. 

Alternatively, one may imagine integrating out the Q(r)-field of the normal region, subject to the boundary con- 



dition Eq. (|15|) and weighted by the action St + Sd- This procedure (we shall describe it in details in section IV) 
leads to the effective action Std of the interface plus diffusive region written in terms of Qs and Ajv only. If all the 
relevant energy scales are less than the Thouless energy of the diffusive metal - the general form such action may take 
is0 

Std = £ li Tt{(Q s An) 1 } , (16) 

i=i 

where 7; are coefficients which depend on the details of the interface (in our case the ratio Gt/Gd)- The largeness 
of the Thouless energy is necessary to disregard the retardation effects and thus possible time non-local coupling 
between Qs and An. Under such condition, the proximity action, Eq. ([l6]), is completely equivalent to those given 
by Eqs. ( p"3| ) and ( |l4| ) upon the proper choice of the set, 7; pd| ]. 



III. THE INSTANTON APPROACH 



We are interested in the limit of strong coupling between the dot and the leads, meaning Gd, Gt > 1 (the weak 
coupling limit may be treated in the spirit of Refs. p2|-^|). For Go 1 the fluctuations of the Q-field around 
its optimal value are suppressed. One may employ therefore the stationary phase treatment of the NLerM for the 
dot-lead interface fl5|j . Taking the variation of the action Eqs. ([l^), ( |l4| ) under the condition Q 2 = 1 one obtains the 
Usadel equation 

2-kvV(DQVQ)-8{v)G t [Q,Qs\ =0. (17) 

This equation is to be solved for a fixed Qs — Qs\&\ given by Eq. (|J) and with the boundary conditions Eq. (|Ts|). 
The solution Q — Q[Q] after being substituted back into the action Eqs. (p~3|) , ( |l4| ) - results in the semiclassical phase 
action S[$\. The later may then be investigated using the instanton approach applied to the $(r)-field. 



A. Zero winding number 

As a warm-up exercise, consider the zero winding number sector of the theory, W a = for a = 1, . . . 2p. Obviously 
it does not produce an oscillatory dependence of F(q), cf. Eq. ([l^), and therefore serves only an axillary purpose. 
The lowest energy configuration in zero winding number sector is &a{i~) — and therefore Qs on the dot is simply 
given by the BCS Ag, Eq. (||). Solution of the Usadel equation ( |l7| ) may be written as 

Q(r) = A N exp {i u(r) 6 <g> a 2 } , (18) 

where it(r) is the normalized "voltage drop" inside the normal region: u(r) = (L — r)/L in ID and u(r) — 
ln(L/r)/ \n(L/d) in 2D. The absolute value of the voltage drop, 9, is coordinate independent diagonal (in replica 
and Matsubara space) matrix. It has a physical meaning of the SC rotation angle of the normal metal in the di- 
rect proximity to the dot. One may substitute the solution back into the action to obtain (for a single replica and 
Matsubara component) 

S a = l [G D 2 - 2G T cos(0„ - 6)] , (19) 


where the subscript "0" stresses that we work with zero winding numbers and 9 n is defined by Eq. (^). For small 
energy, e n A, one has 9 n « n/2 and therefore the corresponding action takes the form 

5 o = [ d2 ~ 2tsin °] 1 (2°) 



G 



where t = Gt/Gd- This action is minimized when 9 = 9(t) satisfies the equation 

9 = tcos9. (21) 

The lowest energy solution of this equation smoothly interpolates between 9 = for t <C 1 and 9 = 7r/2 for t 1. 
Finally, this solution has to be substituted into Eq. (^) to find the action cost, So = So(t), for the zero winding 
number configuration. 



B. Non— zero winding numbers 



To calculate the oscillatory component of the free energy, F(q), one has to consider the <J>-field configurations with 
non-zero winding numbers, cf. Eq. (|l2"|). Consider, thus, the simplest even configuration of winding numbers with 
W\ — 2 and all others W a = in the remaining 2p — 1 replica. With the exponential accuracy it is sufficient to 
consider the "straight" windings: 3>i(t) = 2itTt. The saddle point of the SC Q-field on the dot is given by Eq. (^) 
and takes the form 



/nil/ \ ( fin,m COS 9 m -\ 5, 

Qs {n,m)= , . ft 



,m+2 sin 9 m+ i 
—5 n ,m cos 9 m+ i 



(22) 



where the 2x2 structure is the Nambu space. In all other replica except of a — 1 the Qg-matrix has the form Eq. (|5). 
One may check that the Usadel equation is solved by exactly the same 0{2) rotation as in W — case, Eq. (|l£), 
performed in each of the following 2x2 Nambu blocks 



2u(n+ l,n- 



Qi2(n 
Q22(n 



l,n- 
l,n- 



cos 9 n sin 0.,, 
sin — cos t 



(23) 



the Nambu indexes are explicitly stated on the l.h.s. In this 2x2 block Ajy = 03 and therefore the matrix from the 
r.h.s of Eq. ([23] ) may indeed be rotated into A at as in Eq. (fisl). The only exception are the two lowest Matsubara 
components n — and n = 1, which form 4x4 block of the form 



/On (1,1) Qia(l,-1) \ /costfo sin O \ 

On(2,2) Qia(2,0) cos^ sin6»i 

Q 2 i(-1, 1) Q 2 2(-l, -1) " sin^o -cos^o 

\ Q 2 i(0,2) 022(0,0)/ \ sin^ -cos9i J 



(24) 



The corresponding 4x4 block of the A^r— matrix on the normal lead is the unit matrix. Obviously the unit matrix 
cannot be unitary rotated into the block Eq. (|24|) and therefore solution of the Usadel equation in this block is 
impossible. The difficulty originates from the fact that, due to the random phase of the CB oscillations, the average 
free energy, (F(q)), is not an oscillatory function. We need therefore to consider winding number configuration of 
the form W\ = 2 and W p +i = —2, while all other W a = 0. The contribution to the correlation function from such 
a configuration is proportional to exp{27ri(<7 — q')}, that is the lowest normal harmonic of the correlation function 
(F(q)F(q')). The "dangerous" 4x4 block in the p + 1-st replica is given by 



/Qn(-1,-1) Qi 2 (-l,l) \ /cos0 o 

Qn(0,0) 0i2(0,2) 

32i(l,-l) 022(1,1) sin0 o 

V O2i(2,0) 022(2,2)/ V 



COS 01 



sin Hi 



sin 9<j 

- COS #0 



sin 9\ 
- cos 9\ J 



(25) 



while the corresponding 4x4 block of the A^r-matrix is minus one times the unit matrix. Here as well the unitary 
rotation between the points r = L and r = is impossible. However, if one combines a — 1 and a = p + 1 replica and 
allows rotation between them - then the unitary rotation may be found flEI . Indeed, combining both "dangerous" 
blocks into the single 8x8 block (i.e. combine Eq. (^J) with Eq. ( p5| ) on the dot and unit matrix with the minus unit 
matrix on the normal lead), one readily see that they may be unitary connected (since they possess the same set of 
eigenvalues) . 

The calculations are simplified in the low-temperature case, T <C A, where 9 n w 7r/2. In this case the 8x8 block 
takes the form 



Qs = 



i 













• 


ox I 



(26) 
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on the dot and 



A -=^0lJH0 -lo) (27) 

on the normal lead (oo is the unit matrix in the Nambu space). Here the outer 2x2 unit matrix represents the space 
of n = and n = 1 Matsubara components, while the inner one - the replica space of a = 1 and a = p + 1 (Pauli 
matrices act in the Nambu space) . We seek thus for the solution of the Usadel equation in the 8x8 sub-space, having 
in mind that the solution for all other Matsubara components and replica is exactly the same as in W = sector, 
Eqs. (|l8|)-(pT|). The general solution for r =^ may be written as 

Q(r) = A Ar cxp[z U (r) ( ° f % \\ (28 ) 



where it(r) function was defined after Eq. ( |18[ ) and B is a coordinate-independent Nambu matrix. Employing the 
singular value decomposition, it may be written as 

where U = e iU3 " 3 e iU2 " 2 e iuifT1 and V = e iV3 " 3 e lV2<J2 e™ 1<Tl are SU(2) matrices and 6i < 9 2 are real singular values. 
Substituting the solution back into the action, one obtains (for each of the two involved Matsubara frequencies) 

2S±2 = 1 [Gd(@1 + e!)-G T (sm2u 2 -sm2u2)(cosei-cose 2 )] , (30) 

where the subscript ±2 stays for the corresponding winding numbers and the coefficient two on the l.h.s. reminds 
that two replica were involved. The next step is to minimize the action over 61,2, 1*2 and V2 angles. Three of the four 
equations for the minima have only trivial (parameter independent) solutions: Oi = 0; U2 = — 7r/4 and V2 = 7r/4. 
The action in terms of the single non-trivial angle 02 finally takes the form 

25± 2 = -j- [Bl + 2i(cos 9 2 - 1)] , (31) 
where, as above, t = Gt/Gd- The corresponding saddle point equation for 02 = 62 it) is 

e 2 =<sine 2 . (32) 

For t < 1 the only solution of this equation is 02 = 0, while for t > 1 the angle 82(0 interpolates between zero (for 
t = 1) and 7r (for t > 1). 

C. CB suppression 

We are now on the position to discuss the suppression of the CB. Consider first the component of the correlation 
function (Z p (q)Z p (q')) which is proportional to cos 27r(g — <?'). As was explained above the relevant field configurations 
are those having a single replica with W a = ±2, where a £ [1, p) and a single replica with W a = T2, where a E [p+1, 2p] . 
The corresponding action is given by 45 = 25±2 + (2p — 2)25o; (here So is multiplied by the number of replica with 
zero winding number - 2p — 2, and by factor of two, because two Matsubara components are different between W = ±2 
and W — 0). Taking the replica limit p — > and employing Eqs. (pp|), (|3l]), one finds 

2S(t) = S± 2 - 2S = [Gjftt) - 29 2 (t) + 2<(cose 2 (t) + 2 sin 9(t) - 1)] , (33) 



where 9(t) and 62 (t) are the solutions of Eqs. fl2l|) and (|32| ) correspondingly. As a result, the contribution to the 
correlation function with the unit period has the form: exp{— 45} cos27r(q — q'), where the factor of 4 stays for the 
fact that two replica are involved in the correlation function and the winding number is two. We shall discuss this 



result in more details in section VI, after deriving it using other methods. For the later reference we need to calculate 
t 2 d t (2S/t); employing the saddle point equations ( ^T|) and (H|), we obtain 
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As was explained in section ||, the parity effects sets in at small temperature T < T* « A. Technically it manifests 
itself in the appearance of the odd winding numbers. The calculations for the field configuration with W\ = 1, and 
Wp+i = — 1 are exactly parallel to the one for the even winding numbers. The only difference is that there is the single 
"dangerous" Matsubara frequency, n = 0, consequently there is no outer 2x2 Matsubara structure as in Eqs. ( |2^ ) 
and ( p7| ) (the replica and Nambu structures are exactly the same). As a result, the action is exactly twice smaller 
than for the oscillations with the unit period. We thus obtain for the component of the correlation function with 
doubled gate voltage period: exp{— 2(5 + A(T))} cos ir(q — q'), where the factor of two in the exponent originates 
from the two replica involved. 

Finally, let us discuss the technicalities of the replica limit, p — ► [H| . There are altogether p 2 distinct configurations 
of the winding numbers which contribute to the correlation function (Z p (q)Z p (q')} with the same oscillatory factor, 
say exp{27rz(<7 — q')}. They correspond to the p possible choices for W a = 2 and p independent possible choices of 
W a — —2. There is thus p 2 combinatorial factor in the (Z p (q)Z v (q')) that is cancelled against the same factor in the 
denominator for (F{q)F(q 1 )) = lim p ^ p- 2 {(ZP(q)ZP(q / ) - 1)). 

IV. PROXIMITY ACTION 

In this Section we discuss the derivation and simple applications of the proximity action approach to treat the 
interface between the SC dot and normal lead. This method was recently suggested jll],^5| for the analysis of 
superconductive- normal metal transition in a 2D proximity-coupled arrayj2(|,|27j . Notice that our definition of the 
running parameter £ and charges {7/} differ from those used in Refs. [[ll 25 27J], see Appendix |X| for details. 



A. Derivation of the proximity action 

Our immediate goal is to derive the proximity action, Eq. (|l6|), starting from Eqs. ([l3|) and (|l4|). To this end, let us 
imagine adding the diffusive metal step by step with infinitesimal shells having conductance 6G ^> Gr> > 1. At some 
intermediate step of this procedure one has the tunneling barrier with the part of the normal metal incorporated into 
the set of coefficients 72(C)- Here C = Gd/Gd(0 S [0, 1] is the running parameter such that 7;(0) = Si^Gt/Gd, and 
S(0) = St is the bare tunneling barrier action; while 7;(1) = 7; and 5(1) = Std is the full proximity action of the 
barrier "dressed" with the diffusive metal. At each step one adds another thin shell of the diffusive metal. The action 
of the entire system takes the form 

G 00 hG 
S(C + SQ = — f- 7/(0 Tr{(Q s Q)<} + — Tr{(Q - A N ) 2 } , (35) 

{=1 

where Q is the field right on the boundary between the new shell and the already integrated region. The action of 
the newly added shell has the simple form, (SG/16) Tr(Q — A^) 2 , since SG ^> Go > 1 and therefore the matrix Q 
has to be rather close to A^ (in the leading order in (SG)^ 1 one may disregard all higher order terms (Q — A^) 1 with 
I > 2). The next step is to integrate out the Q field and obtain the new set 7;(C + SQ. We perform this procedure in 
the saddle point approximation. Taking variation of the action, Eq. (|35|), under the condition Q 2 — 1 one finds the 
saddle point (Usadel) equation 

00 

G D WO [(QsQ) 1 - (QQs) 1 ] - SG [QA N - A N Q] = . (36) 
1=1 

In the leading order in Gd/SG <C 1 one obtains for the saddle point solution 

„ 00 

Q = An + 25& ^2 hl [(Qs^n) 1 A n - A N (Q S A N ) 1 ] . (37) 
1=1 

Putting it back into the action, one finally finds 
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OO « * oo 

^ 7/ Tr{(g s A w y} + ^ Yl lk 7nk Tr{(QsA N y +k -(QsA N ) ]l - kl } 

1=1 l.k=l 



(38) 



where we employed the fact that 5G = Gd/$(- One ends up with the following evolution equations for ji(C)' 



^ k=l k=l 



(39) 



To solve this set of equations it is convenient to define the function m(£,8) 
Eq. ( |39] ) takes the form 

1t£ = —uuq , 



£ Z7j(C)sin(/8). With it's help 
i=i 



(40) 



supplemented with the initial condition u(0, 8) = isin8, where, as before, t = Gt/Gd- Employing the method of 
characteristics, one may write the solution of the latter problem in the implicit form: 



«(C,6) = tsin(e-Cu«,9)). 



(41) 



As a result, the function uo(Q) = u(l, 8), that has l-fi as it's Fourier coefficients, may be found as the solution of the 
following differential or algebraic equations 



t l 



d_ 

di 



+ u o d u o = ; «o(6) = t sin(8 - u (Q)) . 



(42) 



The last expression solves the problem of writing the action of the barrier and the diffusive region in the form of 
Eq. (|l(]). Note, that to derive this result we have not used a specific form of Qs and An (other than the fact that 
Q s = K? N = 1). Therefore the action and the coefficients ji are equally applicable to any other setup. The only two 
conditions that were employed are Go S> 1 to disregard localization effects and perform the saddle point calculations 
and i?Th ^> Ec to disregard the frequency term in the action and justify the structure of Eq. ( |l6| ) . The condition 
-^Th 3> Ec also allows to neglect Cooper-channel interaction in the normal conductor, while deriving proximity action 
(cf. |JTT[] for the discussion of general situation). 

For the later use we define 9 = 8(t) = uq(tt/2) and O2 = 02(t) = Uo(t)- According to Eq. (|i^), these two angles 
satisfy 8 = tcos9 and ©2 = t sin 82 correspondingly. Comparing these relations with Eqs. (|2l| ) and (|3^ ) , we conclude 
that the angles 9 and 82 introduced here coincide with the instanton angles introduced in section 



B. Physical observables 

In this subsection we repeat briefly derivation presented originally in |ll[] . 

Normal state conductance. Consider for a moment the dot in the normal state with a small applied bias V(r) = 
V cosuj m T. The dot's Green function takes the form Qs — exp{i$(-r)}Ajv(i" — t') exp{— i$(t')}, where $ = J cLtV{t). 
Substituting such Qs in the action and taking the second variation with respect to V(r) at V = 0, one finds for the 
linear normal state conductance 

00 

G N =G D J2l 2 H ^GodeuviO). (43) 
1=1 

For small 8 one has Mo(8) ~ Q8qUo{0); employing Eq. ([i"2"|), one finds <9eMo(0) = t/(l + t). As expected, one obtains 
the resistance addition rule 

G N ' - G D ' + G~j} ■ (44) 

Andreev conductance. We shall concentrate now on the temperature region A > T > where the SC Q-matrix 
may be written in the form Eq. (^) with 9 n = it/ 2, that corresponds to infinite A (the leading effect of the finite A 
is to renormalize E^} — > E^, 1 + G/2itA plp^|). In this case the Qs matrix is off-diagonal in the Nambu space and 
time-local. Since the Ajv-matrix is diagonal in the Nambu space, only even powers of (QsAjv) may contribute to the 
action. As a result, one finds for the proximity action 
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Std m = ~ E ™ Tr { (Qs An) 21 } 



(45) 



1=1 



£(-l)Sl £ / rfr i • • ■ / dT * e 2i * a(Tl) G(n ~ r 2 )e- 2 ^^G(r 2 - r 3 ) . . . e- 2i *"<™>G(r 2i - n ) , 

i=l a=l n n 



where G(r — r') = — iT / svh-kT(t — t') is the Green function of the normal lead (the Fourier transform of Ajv(e„) to 
the time domain is Ajv(t, t') = G(r — T r )a^). If the external voltage V(t) is applied to the SC dot, the phase is to 
be understood as $ = J drV(r). Calculating the second variation of the above action with respect to V(r), one finds 
for the linear (Andreev) conductance 



G A = G D ^(-1)'(20 2 7 2/ = G D d e u^/2) . 



1=1 



Employing Eq. (W2h for 6 « 7r/2 one finds 



G A 



Gt sin 



(46) 



(47) 



1 + t sin 6 

where 9 = tcos9 = uq(tv/2). One therefore obtains Ga = G^/G D + Q(t) for t < 1 and Ga=G d + 0(1 /t) for t > 1. 



C. Instanton treatment of the proximity action 



To calculate the proximity action on the simplest instanton trajectory < &i(t) = ttTtW and &o,(t) — for a € [2, 2p 
one may employ e.g. the Matsubara basis (see Appendix [X] for calculations in the time domain). Employing Eq 
(modified in the obvious way for arbitrary W), one finds 



(Q s AN) 2l (n,m) = (-iyS n ^ 



sign(e„_^)sign(e r , 






sign(e„ +vv )sign(e„) 



(48) 



in the replica a = 1, while (QsAjv) 2 ' = (— 1) <7o m all other replica. As a result, for all even I — 2k the action in a = 1 
replica is not different from that in the other 2p — 1 replica and therefore does not contribute to the total action in 
the replica limit p — ► 0. On the other hand, for odd I = 2k — 1, there are \W\ Matsubara components, where a = 1 
and a / 1 replica come with the opposite signs. Employing Eq. (p5|), one finds for action on the instanton trajectory 
(in the replica limit) 



Sw(t) = -|W|^E74 fc -2 = ~\W\^- 



k=l 



tt/2 Tr 

rfOMo(O) - J dOu (Q) 

tt/2 



(49) 



where the last equality is a direct consequence of the definition uo(Q) = Y^=i ^li sin(ZO). Employing the differential 
equation (Eg)), one finds 



Gd 
~8~ 



tt/2 tt 

J dQd B ul - J ded e u 2 



7T/2 



ell 



(50) 



Comparing with Eq. (|34|), we find an exact coincidence of the real-space result, Eq. (p3|), and the proximity action, 
Eq. (p)|), 5(t) = 5i(t). Notice that, since the (random) diffusive region was integrated out upon derivation of the 
proximity action, the phase of the CB oscillations is not random. In a sense the proximity action represents a typical 
diffusive region (mesoscopic fluctuations are omitted in the derivation given above). One may therefore calculate 
directly typical Z(q), without resorting to the correlation function. From the two smallest winding numbers W — 1 
and W = 2 one finds Z(q) ss exp{ — (Si + A(T))} coswq + exp{— £2} cos 2irq + . . ., in agreement with the result of the 
previous section. 
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It is worth mentioning here that the general expression for the instanton action in the case of the normal island, 
derived in Ref. fllj], may be presented in the form similar to Eqs. (p3f) and (|4£ 



oo n 

\w\-f- 5> 2fe -i = \w\-2 [ deuo(e) = [e| + 2t( cos e 2 + 1)] 

k=i i 



(51) 



and only even W are allowed for the normal case. It is easy to see, using Eq. (32), that in the entire interval Gd > Gt 
the value of the action S2 is the same as in the tunnelling limit: = Gt/2. As t = Gt/Gd approaches unity, a 
kind of phase transition occurs [^9|, which reveals itself as a non-analytic behavior of (t) at t = 1. 



D. Landauer approach 

Employing Nazarov's techniques one may show that the function uq(&) is directly related to the generating 
function of the transmission coefficients 

N 

\ G D uo(Q) = F(0) sine ee ]T - ^° Sin 2 ° > ( 52 ) 
2 1-T a sm (6/2) 

where T a are eigenvalues of the N x N matrix vt describing transmission between the dot and the lead. Factor 1/2 
on the l.h.s. takes into account that is defined for the two spin components. Performing the Fourier transform 
one finds 

a=l 

where T a — cosh -2 [i a . As a result, the proximity action takes the compact form: 

N . N 

2 



N N , v 

5 T D = -^Trln(l + e - 2 ^Q s A JV ) =--^Trln ^-^(Qs-Ajv) 2 . (54) 



a=l 

In the last expression we have omitted a constant which goes to zero in the replica limit. Exactly the same interface 
action is known in the supersymmetric formalism |3^|| (notice that here we deal with spin 1 /2 electrons) . One may 
calculate now the action on the instanton trajectory to find the renormalized charging energy. Employing Eq. (|53| ) 
and recalling that Si = —(Gn/2) J2k 74fc-2, one finds 

1 N ( T 2 \ 

^--sSK 1 -™- (55) 

This expression is to be compared with the one for the Andreev conductance obtained from Eqs. ( fi"6| ) and (53) [ ]l7| , |l8| ]: 

Q — 1 V ' 

(c.f. with the normal conductance Gn — G]j^2 l l 2 ji — 2^ Q T Q , where 2 stays for spin). Factor 4 on the r.h.s. 
of Eq. ( p6| ) stems from the fact that Andreev particles carry charge 2e and no spin. One may thus identify the 
combination T a = T 2 /(2 — T a ) 2 as the Andreev conductance in channel a. The renormalisation of the charging 
energy is therefore given by 



/v 



= n( 1 - T «) 1/2 ' ( 57 ) 

a=l 

This may be compared with the corresponding factor for the normal spinless particles EJ : J1q(1 — Ta) 1 / 2 . Indeed, in 
our notations the normal instanton action is = (Gu/2) 722-1 = — Y]„ ln(l — T a ); for spinless particles one has 
exp{— 5^/2} for the CB suppression. There is thus a perfect analogy between Andreev and normal spinless charge 
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fluctuations mechanisms, provided that the normal transmission coefficients, T a , are substituted by the Andreev ones, 
T a . 

From Eq. ( |5^ ) one may extract some general results. In the tunnelling limit: T a <C 1 in all channels, one finds: 
Si = Ga/8. In the diffusive metal interface the transmission eigenvalues are distributed according to the Dorokhov 
distribution jl6|,^9|]: P(T) — Gn / — T). The corresponding distribution for the Andreev transmissions is given 
by P(T) = G_d/(8T-v / 1 — T), as a result Ga = Gd- For the typical action of the diffusive interface one finds therefore 
(Si) = n 2 G A /32. 



E. Fluctuation determinant and summation over many-instanton configurations 

We now expand the action, Eq. (fl5|), to the second order in deviations from the instanton trajectory = 
wTrWa + (5$ a (r), where <5<I> a (0) = <5<E> a (,0). After diagonalization of the resulting quadratic form one finds the 
following spectrum of the small fluctuations O,0,m[: 



TT 2 T 

"Be 7 ' 



^(\n-W a \ 



W a \-2\W a \) , 



(58) 



where n G [—00,00]. For T <C Ec there are 2|T / F a | + 1 (almost) zero modes labelled by n < \W a \- One of them 
n = is the trivial shift of Qa(T~) by a constant, whereas the remaining 2|W / a | zero modes are associated with the 
deformation of the instantons. The general solution of the saddle point equations (for T <C Ec) known as the 
Korshunov instanton |32| may be written as 



\w a \ 

n 

fe=i 



1 



sign(VK a 



(59) 



where Zk is a set of \W a \ complex numbers (instanton coordinates) which parametrize the 2|Wa| dimensional zero- 
mode manifold. The instanton coordinates, z^, are complex numbers from inside of the unit circle. One may thus 
characterize any deviation <5 < I , q (t) by 2| W a \ zero-mode coordinates, z^, and remaining transversal modes having non- 
zero masses, Eq. (^8|). The Jacobian of such transformation may be found in the standard way j3^| and is given 
by 



J 



1 



\w a \\ 



det 



1 



1 



z l z k 



(60) 



where I, k € [1, |W |]. The fact that the Jacobian vanishes if two of the coordinates coincide signals the repulsive 
interaction between the instantons. 

The Gaussian integration over the massive degrees of freedom normalized by those in the zero winding number 
sector results in 



n (Ai 0) A) 
n=i 

n (a™/*) 

n=\W a \ + l 



Ga 
2tt 



\W a \ 



cxp UW a \ln 



( G A E C \ 
\ 2ir 2 T J 



(61) 



up to the terms of the order T/(GaEc) -C 1- Combining together all the factors, one finds for the one instanton, 
W a — 1, trajectory 



Zi 



d 2 z 



1 



Ga 
2tt 



-{Si 



An(G A E c /2^T)} _ g 



^£ In — 



2T 



T 



(62) 



The cutoff of the logarithmically divergent ^-integral originates from the following consideration. In presence of the 
charging energy term the Korshunov instantons, Eq. (|5^), are not true zero-modes of the action. The z-dependence 



of the action brings the factor exp{— (T/2Ec)\z\ 2 /(l 



'■)} to the integral. As a result the effective integration 



range shrinks to \z\ < 1 — T/Ec < 1. The renormalized charging energy Ec is defined as 



E ( 



E C G 2 



2tt 2 



A e - Sl 



(63) 
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It is also instructive to look at the W a — 2 contribution to the partition function 



^2 = ^iTtq I 



2! 



2! 



bd 2 z\d 2 Z2 



2T 



hi 



(i-N?)(i-N!) {i-z lZ * 2 ){i-ztz 2 ) 



2 #C 7T 2 



^ \ p -{S2-2\n(G A Ec /2tt 2 T)} 

2^ ' 



T 6 



(64) 



The leading power of the large logarithm comes from the diagonal term of the Jacobian, Eq. (|6C|). This term 
corresponds to the non-interacting instantons approximation, which is justified therefore by the large parameter 
ln(Ec/T) ^> 1. Notice, that the (repulsive) interaction correction (the off-diagonal part of the Jacobian, Eq. (|6C|)) 
comes with In and not In 1 , as may be naively expected. 

Due to the large logarithmic factor, originating from the integration over sizes of each individual instanton (i.e. 
neglecting their interactions), one may treat positive and negative instantons (anti-instantons) on the equal footing. 
Strictly speaking, a combination of instantons and anti-instantons is not a saddle point solution. However, relative 
weakness of instanton interactions (no log-factor in the second term in Eq. (64)) makes it possible to consider an ideal 
gas of instantons and anti-instantons. The contribution to the partition function with a given winding number W is 
given by all configurations having m + \W\ instantons and m anti-instantons in an arbitrary order. Taking only the 
terms with the leading power of the logarithms (diagonal term of the Jacobian) one finds Q 



Z(q) = Z Q 

W=-o 



1 



oo 

-KiWq \ 

^ m!(m+ \W\)\ 

rn—0 v ii/ 



in 



Ec 
T 



2m+\W\ 



(65) 



This ugly looking series may be summed up into an unexpectedly simple expression: 



Z(q) = Z exp | In ( ) o >M M \ 



(66) 



The simplest way to check it is to expand back Eq. (|36j) in a double series in powers of Eq and e l7!q . The role of the 
anti-instantons, therefore, is to convert cos(Wirq), which may be expected from the instantons only, into cos w (irq) . 

Under the condition Ec < T <C Ec the interaction between instantons can be treated perturbatively; for the 
lowest-order correction we use the last term in Eq. ( |64]) to find for the gate voltage - dependent free energy 



F(q) = -E c 



111 



Ec 



cos(vrg) - ^^-cos(27rg) 



(67) 



We retain the second term in (|6?| ) since it will be important at temperatures comparable with the parity-effect 
temperature T*, see Sec. VI. It also demonstrates that the instanton-instanton interaction corrections are small only 
if E c < T. 



V. RG TREATMENT OF THE QUANTUM FLUCTUATIONS 



We next approach the problem of the CB from the different perspective. Instead of calculating the action on 
the instanton field configurations, we shall perform the RG analysis of the coefficients "fi upon integrating out fast 
fluctuations of ^(t). To this end we write ^(t) = $ s (r) + ^(r), where superscripts s and / stay for slow (ui m < 57) 
and fast [u) m > fi) component of the field correspondingly. The running cutoff fl runs from ~ GaEc down to the 
lowest Matsubara frequency £1 ~ T. We next substitute the $-ficld into the SC proximity action, Eq. (45), expand 



to the second order in $^ and integrate out the fast field fluctuations. The bare propagator of the -fields may be 



read out from the action Eq. (45) 



($f(u m )$l(-w m )) 



\vm\G A ((n) ' 



(68) 



where Ga((q) = Gd X^i(~l)'(20 2 72z(Cfi) ls runnm g value of the Andreev conductance, and £n = ln(G aEc) 



(Note that this definition of C differs from that used to derive the proximity action in Section IV.) The further 
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calculation is essentially similar to that in the Keldysh formalism [E7| . One obtains the following set of the evolution 
equations for ji(Cn)- 

where we have omitted the subscript f2. Since only the even coefficients are involved it is convenient to define the 
function u(£, 9) = [u((, O) — tt — 0)]/2 which contains only even harmonics. In terms of this function the RG 



equations (69) take the form 



uc = -^j[fi(e)taa0] e , (70) 

where the Andreev conductance is Ga(C) = GdUq(C,tt/2). The initial condition for Eq. fl7C| ) is u(0, 0) = Uo(0), 
with the uo(0) function evaluated above. Let us mention for completeness, that in the normal case the corresponding 
RG equation takes the form = (2/Gj\r(£))[u(0) cot(0/2)]e, where the running normal conductance according to 
Eq. © is given by G N (C) = G D ££ x ^7/(0 = G D u e (C, 0). 

The key parameter that determines the strength of phase fluctuations is the Andreev conductance Ga(Q- We con- 
sider the large bare value, Ga(0) S> 1. Upon integration over fluctuations of $(r), the effective value of Ga{Q decreases 
(the physical reason of this phenomena is a partial loss of the phase coherence between multiple Andreev reflections), 
and becomes comparable to unity at some time scale il^ 1 = (GaEc) -1 ^" ■ At longer times coupling between the 
dot and the lead is weak, thus the phase fluctuations grow and autocorrelation function C(t) = (e l (*( T )~*(°))) decays 
fast. It is natural to associate f2 c with an effective Coulomb energy of the island connected to the wire, Ec- We will 
show now that the estimate for Q c that follows from the RG equation coincides (within exponential accuracy) with 
the results of the instanton analysis. 

Equation ( |7C| ) may be solved with the method of characteristics. We define the function u(0) = u(Q) tan© and 
then introduce new "space" and "time" variables: £ = — In sin© and t which is defined via the relation 

dQ^dr. (71) 



Ga(0 



In terms of these new variables the RG equation (f70| ) takes the simple form: v T = , yielding v(£, t) = «o(£ + r). As 
a result, 

P T POS (~~) 

u(C(r),0) = — u (arcsin(e- T sin©)) , (72) 

VI - e~ 2T sin 2 



The dependence £(t) can be now found from Eq. (|7l| ) with the Andreev conductance given by Ga(C) — Gduq((, tt/2). 

Our goal is to find a scale £ c , where the Andreev conductance Ga becomes small. To find £ c we notice that in terms 
of the RG time t the smallness of Ga means the limit t — *■ oo. Then, using definition of r given by Eq. (|7l|), we find 



oc 



tt/2 

f u (arcsine- r ) G D f 
Cc = -— Jdr - ~— J dO uo(0) . (73) 



Recalling the definition of Uo(0) and comparing this result with Eq. ( fl9| ) we conclude that ( c = Si. As a result, the 
amplitude of the lowest oscillatory component (~ cos nq) of the free energy, from the RG point of view, is given by 
the precisely the same exponential factor as in the framework of the instanton calculation. 

Yet another instructive way to approach the problem is to follow the renormalization of the instanton action. 
According to Eq. ©, 5^(0 = -\W\(G D /2) Efc74fc- 2 (C) - -\W\[G D /A) ^^dQu^Q). Employing Eq. ©, one 
finds 

§^^("(0)tanO) „=-\W\, (74) 
«C G A (Q e-^f 

where in the last equality we have used that Ga[C) = Gduq((, 7r/2). (In the normal case one finds dS^/d( — — \W\/2; 
only even W are allowed.) From the dimensional analysis one expects F(q) ~ il exp{— S±} cos nq + . . ., where 
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fi ~ EcGa is the cutoff energy. We may now renormalize down the cutoff energy fi — ► fi(C) = fie simultaneously 
following renormalization of the action, S\ — > <Si(C)- According to Eq. (f74|), <S>i(C) = <Si — C- As a result, 

0(0 e~ 5l(0 = const. (75) 

This means that the lowest oscillatory component of the free energy is expected to be temperature independent, in 
agreement with Eq. (|67j) (up to log-factor, originating from the zero modes). The instanton calculation is done at 
C = 0, alternatively the RG calculation presented above looks for £ c , such that Sw(Cc) ~ (the scale where the 
instantons do not cost anything). In view of Eq. ( f75| ) it is not surprising that they yield the same result. Moreover, 
one can renormalize down to any intermediate scale C (such that Ga(() > 1) and then calculate the instanton action 
with the current parameters {ji(()} ~ the result is still guaranteed to be the same. Notice also that the renormalized 
action Sw(C) — S\y — \W\ \u(GaEc/T) coincides with the result that comes from the instanton action together with 
the massive Gaussian fluctuations, c.f. Eqs. (|6l|)-(|64|). The effect of the zero modes (In Ec/T in the prefactor) is 
missed in the one-loop RG calculation. 

Finally we comment upon the comparison of the approach used in this Section with Ref. |p7f where the same 
kind of RG approach was employed in the case of 2-dimensional normal conductor with very low Thouless energy, 
-&rh *C Ec- In the last case the RG procedures involves simultaneous integration over phase fluctuations ^(t) and 
diffuson/Cooperon modes in the normal conductor. As a result, the full functional RG equation contains three terms 
in the r.h.s.: one coming from Eq. (|io|), the second is similar to Eq. (f70|), and the third contribution accounts for the 
Cooper-channel repulsion in the normal conductor (cf. ]ll]^5| ) . The resultant RG equation cannot be reduced to the 
set of even harmonics u only; we are not aware of any method to solve it apart from "brut-force" numerical integration. 
Considerable simplification of the problem treated here stems from the fact that, due to the condition i?Th ^ Eq 
all diffuson/Cooperon modes may be integrated out before (means at larger energies than) the phase fluctuations 
$(t) become relevant. This procedure leads to the action functional Std that depends upon <&(t) trajectory only, 
as defined in (|45|) . Thus the results presented in this Section can be considered as "minimal generalization" of those 
obtained in Ref. ||. 



VI. DISCUSSION OF THE RESULTS 



The main result of our study, given by Eq. (|67j) , applies in the temperature range 

E c <T«min{E c ,T*}, 



(76) 



where T* is the parity-effect temperature given by Eq. (|lO|). The condition T <C T* ensures that both even and 
odd winding numbers contribute to the partition function, since A(T) < 1, c.f. Eq. (|). The two lowest oscillatory 
components of the free energy are given then by Eq. (|67j ) . At higher temperature T ~ T* an addition of odd number of 
electrons to the dot becomes possible, changing the relative amplitude of the harmonics. On the level of the partition 
function the components with odd winding numbers acquire the temperature dependent factor exp{— A(T)}. As a 
result, one finds for the gate voltage dependent free energy: 



F(q) = -E c 



• MT > 1" ( ^f) coBfo) + ((1 - e-W)ln 2 (^f) - n - 



Ec 
4T 



cos(27rq) 



(77) 



Finally at T > T* the parity effect disappears and only the normal (unit) oscillation period remains in the free energy 



no) 



-E c - 



ec in 2 m i 



AT 



cos(2irq) 



(78) 



Notice that if In Ec/T* > 7r/\/6 there is the sign change of the cos(27rq) component at T w T* . Moreover, unlike 
the low temperature case, where the amplitude of the dominant harmonics was weakly (logarithmically) temperature 
dependent, Eq. (|67|), at larger temperature there is the stronger dependence (~ T _1 In 2 (Ec/T)) of the amplitude. 
The characteristic crossover temperature is determined from the relation A(T^) « \n.(T* / Ec) and is given by 



A 



In 



(79) 
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FIG. 2. (a) Normalized actions 8Si/Ga (full line) and 2S2 /Gn (dashed line) as functions of t = Gt/Gd- Both graphs 
interpolate between 1 in the tunnelling limit and 7r 2 /4 in the diffusive limit, (b) Andreev (full line) and normal (dashed line) 
conductances normalized by Gd as functions of t. 



where we have used the fact that Si ~ Ga- 

We emphasize an interesting feature of the weak Coulomb blockade compared to the usual one: although the role 
of effective Coulomb energy is taken by the effective energy Ec <C Ec, the oscillation amplitude depends relatively 
weakly upon temperature within the range given by Eq. ([76]). This is to be contrasted with the usual case of an 
island connected by highly resistive contacts: at temperatures T > Ec oscillation amplitude vanishes exponentially 
fast with T/4E C - 

Next we shall discuss the quantitative value of the renormalized charging energy Ec, which determines both the 
magnitude of the CB and the region of applicability of our results. To this end we need to evaluate Si(t) (where 
t = Gt/Gd), which is given by one of the two equivalent expressions, Eqs. ([33]) or (f49|). In the two limiting cases 
t -C 1 (the tunneling barrier limit) and t > 1 (the diffusive metal limit) the answer may be found analytically. One 
obtains 

2tt 2 Ec Ga ( 1 ; in the tunneling limit, . , 

G\E C = ~~8~ { tt 2 / 4 ; in the diffusive limit, > 

where the Andreev conductance is given by Eq. (|47|). Recall that for a normal dot the corresponding exponent is 
given by Gat/2 in the tunnelling limit and Gat7t 2 /8 in the diffusive limit. We observe that upon the same dissipative 
conductance the SC dot exhibits factor of four (in the two limiting cases) smaller action. Therefore one may observe 
a sizeable CB in the superconductive state, while in the analogous normal dot the CB is practically suppressed. The 
crossover behavior for the SC action Si (t) , the normal action S™ (t) and the Andreev conductance Ga (t) as function 
of t = Gt I Gd is depicted in Fig. ||. Dependence of the same quantities on the resistance G^, 1 of the normal region 
at fixed Gt is demonstrated in Fig. [|. 

In the case Ec "C the crossover between T < and T > T* regimes leads to the sharp drop of the oscillation 
amplitude by the factor In -pf <C 1 The same drop in the residual Coulomb blockade may occur as function of 
magnetic field, due to the suppression of the A and thus of by the magnetic field. This may result in a strong 
negative magnetoresistance at low temperatures of a granular media made out of small superconductive grains (cf. 
e.g. p4| for the discussion of relevant experiments). This effect can occur if the low-temperature oscillation amplitude 
Ec is sufficiently large to destroy Josephson coupling between grains, so transport of Cooper pairs between grains is 
blocked by Ec- However, quantitative theory of such an effect is still to be developed. 

Another open question concerns the behavior of F(q) in the zero temperature limit, T < Ec- As temperature 
decreases below Ec, all approximate methods, we used, run out of their applicability range: the renormalized conduc- 



17 




G T /G D 

FIG. 3. The same as Fig. |2|but normalized by Gt- Notice that the normal action exhibits a continuous "phase transition" at 
Gd = Gt- The Andreev action has a maximum at Gd ~ 0.8Gt which is a direct consequence of nonmonotonous dependence 
of the Andreev conductance on the resistance of the diffusive metal. 



tance G a drops below unity, fluctuations become strong and there is no apriori reason to treat them within Gaussian 
approximation. Whereas the overall scale of the oscillation amplitude is probably given by its first harmonic, Eq. ( |67| ) 
at T ~ Ec, and is of the order Ec m t^, the precise shape of the oscillations is still to be determined. It is possible 

Ea 

that in the T — > limit the function F(q) becomes nonanalytic at the degeneracy point q — 1/2. An extreme case 
of such a nonanalytic behavior - finite steps in dF/dq at half-integer q - was found in Ref. H] where the Andreev 
conductance was completely neglected. 
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APPENDIX A: INSTANTONS IN THE IMAGINARY-TIME PROXIMITY ACTION 



1. Transition to imaginary time 

The method of the multicharge proximity action was initially developed Jll],|25| in the Keldysh real-time represen- 
tation. Here we show how the Matsubara proximity action, Eq. (|l6|), may be obtained from its Keldysh analog by 
analytic continuation to imaginary time. Since this Appendix serves illustrative purposes we will consider the simplest 
case T = 0. As discussed in the bulk of the paper, the form of F(q) is unknown at T = since at T < Ec the dilute 
instanton gas approximation fails and one has to consider an interacting instanton liquid. Nevertheless, the overall 
scale of oscillations in F(q) may be inferred from the single instanton action which is temperature independent and 
may be calculated at T = 0. 
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FIG. 4. The initial Keldysh contour (a) and its deformation to imaginary times (b). 



To proceed we need to establish a correspondence between our notations and those adopted in Refs. |35| , |i l|j25[ . 
The latter notations will be designated by a prime. Firstly, our conductance quantum Gq = e 2 /2irh is different from 
G'q = e 2 /% used in those papers. Secondly, the running parameter £ and charges {jn} are related by 



C In G D 



C In 87rV' 

Thus, the Keldysh multicharge proximity action may be written in our notations as |3£ 

.G 



SUc [<!>(*)] = £ 7n(0 Tr(Q s A w )" , 



(Al) 



(A2) 



where the trace is taken over time indexes, Nambu and Keldysh spaces. Here Q$ is a matrix in the superconducting 
dot acting in the Nambu space as 



e 2l *C) 

e -3i»(t) 



(A3) 



where $ = diag($>, <&<) is a matrix in the Keldysh space, with <fr>(i) and &<{t) being the fields residing on the 
forward and backward branches of the Keldysh contour. In the normal lead, An(E) = Aq(E) <g) 03, where 03 is the 
Pauli matrix in the Nambu space, and the matrix An(E) acts in the Keldysh space: 



A (E) 



1-2/ -2f 
2(/-l) 2/-1 



(A4) 



f(E) being the distribution function. 



Tracing over the Nambu space reduces the action (A2) to the form 
. Gd 



n=l 



]T(-l) n 72„ / dh... dt 2n tr K e^^AoCii - t 2 )e- 2 ^A (t 2 - i 8 ) . . . e - 2l * ( * 2 " ) A (t 2 „ - h), 



(A5) 



where the trace is taken only over the Keldysh space. 

Transition to the imaginary time is achieved by the deformation of the Keldysh contour C which initially run over 
the time axis from —00 to 00 and then in the backward direction. The desired deformation introduces a vertical 
segment of length (3 at some time t = to so that the contour originates at —00 + i/3/2, runs through the points 
to + i/3/2, t to 00 and then back through to, t a — i/3/2 to — 00 — ifi/2, see Fig. |j. For the purpose of the evaluation of 
the instanton action the choice of the point to is somewhat arbitrary. Note, however, that were we to consider, e.g., 
the correlation function of the form (N(q)N(q')), it would be necessary to allow for adiabatic evolution of the gate 
voltage q = q(t), q' = q(t') and to introduce two vertical segments at times t and t'. 
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The analysis is simplified in the zero-temperature case. Choosing the position of the vertical segment at to = 0, 
we deform the forward segment (—00, 0) of the Keldysh contour to the upper imaginary half-axis and the backward 
segment (0, —00) to the lower imaginary half-axis. The remaining horizontal appendix of the Keldysh contour running 
from to 00 and then back to t = can be neglected as it does not contribute to thermodynamic quantities. 

In order to continue the function Ao(t) to the complex plane in time one should employ its analytic properties [ B7[ . 
In deforming the contour to the imaginary axis the element AJ 2 (r) enters only with r > 0, and the element Aq^t) 
enters only with r < 0. Under this condition the function Aq(t) can be substituted by [^8| 

A (r) - F(r) ( \ \ ) , F(r) = --L. (A6) 
Eq. (|A6|) solves the problem of analytic continuation of the Green function Aq (t) . 



The trivial matrix structure of Eq. (A6) ensures that each term in Eq. (A5) can be written as a multiple integral 
of the single function 

= ( HI' f rT> n' (A7) 
v ' \ <J><(t), for r < 0, y ' 

defined on the whole imaginary axis. The infinitesimal element dt is transformed to 

dt-l'J*"' f rT> °' (A8) 

[idr, for r < 0. v ' 

Finally, we obtain ]36| ] 

00 r.00 

Vox[$W] - ~ 72n / dn... dr 2n e «*(n)-2i*(T i i)+...-«»(^.) p{n „ T2)f{t2 _ ^ _ _ T ^ (Ag) 
that coincides with Eq. (|4q) in the zero-temperature limit. 



2. Instantons in time domain 



Here we calculate the action on the trajectories for which e 2l *( T ) is an analytic function in the upper half-plane. 
For such a solution, the instanton's winding number irW = A$ = J dT(d$>/dr) is positive. Integration over t%, £3, 
. . . , t 2n -\ is performed with the help of 



J-' 



2i$(n)_ 



TO -TlTl- T 2 



2i$(r ) _ 2i*(r 2 ) 

= -i- 7re 2 ^5(r -T 2 ), 

TO - T 2 



(A10) 



that is obtained by regularizing 1/t = [l/(t — iO) + l/(t + i0)]/2 and making use of analyticity of e 2i *( T ) in the upper 
half-plane. As a result, the integral in Eq. (|A9|) is transformed to 



n 

fc=i 



dm. 



2fc e -2i$(T 2fc ) 



n 

fc=i 



r 2(fe-l) — T 2(fe+1) 



7re 



2l * (T2fe) <5(T2( fe -i)-r 2(fe+1) ) 



(All) 



where To = r^n and T 2 ( n+ i) = r%. Expanding the last product and combining similar terms we rewrite it as 

(-1)"5(0) + ]T C™(-l) n - m K m , (A12) 

m— 1 

where 



/OO "t / , 



2fe e -2i*(r 2fc ) 



fc =1 T 2(fc-1) - T2(fc+1) 



(A13) 



The integrals for K m are calculated recursively with the help of the relation 
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J 2i$(n) _ „2i$(r 2 ) „2i*(r 2 ) _ „2i$(t 3 ) p 2i*(n) _ „2i$(r 3 ) 

£^ e -2»(n)f £ e - = 2i- . (A14) 

TT T\—T% t 2 — r 3 n — T3 

Thereby we get l£" m = 2 m W. F or anti-instantons with IT < the same analysis yields K m — 2 m \W\. 
Now summation in Eq. ( |Al2j ) becomes trivial and we obtain for the instanton action: 



smt)] = -- f £ i2n{(-irm + \w\[i - (-1)-]}. (Ms) 

n=l 



Eq. (A15) formally contains the divergent <5-function of zero argument. This part of the answer gives the action of 
the non-instanton configuration $ = $0 = const and thus drops from the difference Sw = S[$(t)] — S[$o] which 
coincides with Eq. (E9J) obtained in the frequency domain. 
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